UNIFORM BOUNDS FOR EIGENFUNCTIONS OF THE 
LAPLACIAN ON MANIFOLDS WITH BOUNDARY 
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Abstract. Let u be an eigenfunction of the Laplacian on a compact mani- 
fold with boundary, with Dirichlet or Neumann boundary conditions, and let 
— be the corresponding eigenvalue. We consider the problem of estimat- 
ing maxji// u in terms of A, for large A, assuming Jj^j v? = 1. We prove that 
maxjv/ li < Cjif-^'-""^^''^! which is optimal for some M. Our proof simplifies 
some of the arguments used before for such problems. We review the 'wave 
equation method' and discuss some special cases which may be handled by 
more direct methods. 



Let (M, g) be a smooth, compact Riemannian manifold of dimension > 2, with 
smooth boundary dM. Let A denote the Laplace-Beltrami operator on functions 
on M . Consider a solution of the eigenvalue problem, with Dirichlet or Neumann 
boundary conditions. 



The subscript p will always indicate the LP(M) norm. In this paper we consider 
the problem of bounding ||u||oo = inaxa;gM |u(a;)| in terms of A, for large A. The 
size of ||m||oo may be considered as a rough measure for how unevenly the function 
u is distributed over M: If u is 'small' (say o(l) as A ^ oo) on a 'large' set S (say, 
of area \S\ = \M\ — 5) then the condition (||) forces ||m||oo to be at least of order 
therefore upper bounds on ||m||oo imply lower bounds on 5, i.e. the area of 
the set where u is 'concentrated'. It may happen that 5 ^ as A ^ oo, for certain 
(sequences of) eigenfunctions, for example on the sphere or the disk. Also, upper 
bounds on ||m||oo yield upper bounds for the multiplicities of the eigenvalues of A. 

In addition to single eigenfunctions we also consider sums of eigenfunctions of 
the form 



for finite intervals / C M. Here Ui, U2, . . . is any orthonormal basis of eigenfunctions 
with eigenvalues Ai < A2 < . . . ^ 00. (The sum is independent of the choice of 
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1. Introduction 





basis.) A simple argument using a suitable version of Sobolev's embedding theorem 



(see, for example, [H6r85|, Thm. 17.5.3) shows that 



(3) supu[o,A](x) <CA"/2. 

X 

(Here and everywhere C will denote some constant only depending on (M, g).) 
Together with Weyl's law 

(4) #{j : A, < A} ^ 7mA" as A ^ oo 

(with = (27r)-"voleuci(-B")vol(M), B"=unit baU in R") this shows that for 
each X the average size of ui{x),U2{x), ... is of order 0(1). However, if M is a 
sphere or a ball, for example, then there is a subsequence Uji of eigenfunctions with 



> cA^" for some constant c > (see Section 2.3). Our main result 



shows that this is the worst possible rate of growth: 

Theorem 1. Let M" he a compact Riemannian manifold with boundary. There is 
a constant C = C{M) such that any solution of (0) satisfies 

(5) llulloo < CA("-i)/2. 

From the theorem one easily derives the following corollary. This is well-known 
(see 'Related Results' below). 

Corollary 2. Under the same assumptions as in the theorem, the multiplicity of 
the eigenvalue A^ o/— A is at most CA"^^. 

The theorem is already non-trivial in the case of a domain in with the Eu- 
clidean metric: 'Interior' estimates are simpler in this case due to the constant 
coefficients of A (see Section p^ ), but the boundary must still be curved, and this 
causes the main difficulty. While the bound (||) is optimal for balls and spheres it is 



far from optimal for a torus and a rectangle (see Section 2.2), and possibly for more 



general M under assumptions on the curvature. For example, in the case of negative 



curvature, the results of Berard p3er77 imply a bound of CA'"~^''/^/log A. The 
first time that the exponent in (^ was improved in any case of nonzero curvature 
was in the paper ( pS95| ]) by Iwaniec and Sarnak, for certain arithmetic surfaces M. 
In the opposite direction, one may ask for which M the ©(A*^""^-*/^) estimate may 



not be replaced by o(A'" ^)/^). See |SZ01| for results in the case without boundary. 



We do not analyze the dependence of the constant C in (|^) on the metric. The 
methods used imply that it is bounded in terms of a finite number (depending on 
n) of derivatives of the metric and of the geodesic curvature of dM. 

Outline of the proof of Theorem The idea is to use the standard wave kernel 
method outside a boundary layer of width A~^ and a maximum principle argument 
inside that layer. 



Let us first recall the wave kernel method (cf. |H6r68, Sog93|). Certain weighted 
sums over many eigenfunctions turn out to be easier to estimate than single eigen- 
functions, since they have a local character. More precisely, given e > 0, choose a 
Schwartz function p such that 

P>0, P|[o,i] > 1, supp(/5) C (-e,e). 

Here p{t) — e^**'*'p(A) dX denotes the Fourier transform. Proving the existence 
of such a p is an easy exercise. 



If ui, U2, ■ • • is an orthonormal basis of real eigenfunctions of the Dirichlet Lapla- 
cian, with eigenvalues Ai < A2 < . . . — s- oo, then we consider the sum, convergent 
by i), 

(6) ^p(A- Aj>j(x)Mj(y). 

3 

Actually, this is the integral kernel of the operator p(A — V— A), but we won't use 
this fact. If we write, using Fourier's inversion formula. 



1 r°° 



= [p(i)e-*^^r(A) 

= 2[/5(t)cos(tA,)]-(A) -p(A + A,), 

(the superscript " denotes inverse Fourier transform t ^ \) then we get the impor- 
tant identity 

(7) - ^j>A^>3 (y) = mt)K{t, X, y)] - (A) + 0(A— ) 
j 

where the error term is small by and rapid decay of p, and 

(8) K{t,x,y) = KM{t,x,y) = Ycos{tXj)uj{x)uj{y) 

j 

is the wave kernel, i.e. for each fixed y ^ M (the interior of M) it is the solution of 

(— - A,)K{t, x,y)^0 in Rt x M, 
K{0,x,y) = Sy{x) 

(9) ^^K{t,x,y)^Q 

K\x(zdM = 0. 

The convergence of (^ is in the sense of distributions (i.e. weakly) in t, for each 
fixed as the argument leading to (^ shows, for example. 

The point of (^ is that K may be analyzed directly from (^. In particular, 
solutions of the wave equation have finite propagation speed (see [II6r85|, Lemma 
17.5.12, for example, for the easy proof using energy estimates); for K this means 

suppiCc {{t,x,y) : dist(a;,?;) < \t\} 

and that K{t, x, y) depends only on the data (i.e. M and the metric g) in Bt{x, y) := 
{z : dist {x,z) + dist {y,z) < \t\}. Therefore, (Q) shows that the sum (||) depends 
(up to an error 0(A^°°)) only on the data in B^(x, y), and this is the local character 
mentioned before. 

We will use (^ only on the diagonal, i.e. for x = y. From the assumptions on p 
we get 

(10) U[x-i,x]{xf <2[p{t)K{t,x,x)Y{\)+0{\-n. 
so the theorem would follow from an estimate 

(11) I j e''^p{t)K{t,x,x)dt\ < CA"-i. 
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Since p is smooth, this is a statement about the singularities of iiT as a distribution 
in t, uniformly in the parameter x. 

As a first step one may now obtain the interior estimate 

(12) W[A-i,A] [x] < C,A("-i)/^ if dist {x, dM) > e. 

In the case of a domain M C M" one has i^M = K^n for |t| < 2 dist {x,dM) by 
locality (finite propagation speed), so (|T^ ) follows easily from ([l0| ) by using the 
explicit expression for the Euclidean wave kernel in terms of an a;-space Fourier 
transform, see Section 2.4. In the case of a nonflat metric, K can be analyzed near 
t = Q hy the geometric optics approximation (see [Hor68|,|H6r85|) or the scaling 



technique introduced by Melrose in |Mel84|, and this gives the interior estimate in 
that case. An other m e thod to obtain (|12|) uses the Hadamard parametrix for the 
resolvent, see [Ava5C], | Sog88| , for example. 

Since |u(a;)| < u^\~i x\{x) for solutions u of (|l|), ( |l2| ) gives in particular 

(H) \u{x)\ < C,A("-i)/2 if X e M,. 

The claim of Theorem | is that in (|l|') may be chosen independent of e. 
When X approaches the boundary in the argument above, one has to either shrink 
the support of p - which means increasing its maximum since J p — p(0) = 1 
(and this gives only Cc < Ce~^^^ in (p^)) - or analyze the wave kernel in the 
presence of boundary conditions. The difficulty with the latter is that the geometric 
optics construction becomes much more complicated near the boundary because of 
diffraction and multiple reflection of geodesies, and a satisfactory parametrix for the 
wave equation has so far only been constructed near points where the boundary is 
strictly convex or concave (see [MT|, [Mcl80 ). This parametrix was used in ]Gri92| 
to prove Theorem ^ near concave boundary points (e.g. near the inner circle of 
an annular domain). On the diagonal x — y the singularities of K have been 
analyzed in suffic ient d e tail near an a rbitrary smooth boundary by Ivrii, Melrose 
and Hormander ([ [vr80 Mel84 , HorSSj) to allow us to obtain the desired estimates 
at points x outside a boundary layer of width A~^. In the boundary layer, a simple 
maximum principle argument then completes the proof, because there u can have 
at most one oscillation in the direction perpendicular to the boundary. 

Related results. Estimates of ||u||oo for large A are closely related to asymptotic 
improvements over the bound (^) on uy^^ x^ (which is often referred to as the 'spectral 
function of A on the diagonal'). That improvements might be possible is suggested 



by the observation that 



M 



#{j : \j e /}, and Weyl's law 



Carleman 



Car35| was the first to prove the interior pointwise asymptotics corresponding to 
(||) (as A oo) 

(13) U[o,A](2^) =7'A"/^+o,(A"/2)^ for dist (2;, aM) > e 

(for domains in M"; for manifolds see |MP49(| , |Gar53| ). Here 7' = (27r)-"/ Vvoleuci(B"). 
The error term was improved to 

(14) ii[o,A](a;) =7'A"/2 + Oe(A(""^^/^), for dist (a;, 9M) > e 



(|Ava52, Lev53 | in R" , |Ava5(;] for manifolds, |H6r68| for higher order operators on 



manifolds; see ^ ^ 

established by writing v? 



SV96 | for further improvements). The connection to Theorem^ is 



[A-i A] ~ ""[0 A] ~ ""[0 A-ij- Then ( [14D gives immediately the 



interior estimate (|12[), and also shows the optimality of the power (n — l)/2 in (|12[), 
for any M (as opposed to (|2[) which is optimal only for some M). 
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It follows from the very precise and general results in Ivrii's book (| Ivr98| ) that 
the e-dependence in (|lj) may be removed, and hence that our theorem even holds 
with u replaced by U[a-i,a]- Our main point here is the simplification of the ar- 
guments near the boundary. However, not all difficulties can be avoided: We still 
have to refer to results on the boundary parametrix in the treatment of the layer 
{x : A-i < dist {x,dM) < 1}. 

If one considers ||u||p for p G (0, oo) instead ofp — oo then one obtains interesting 
phenomena related to the 'restriction theorem for the Fourier transform' of Stein- 



Tomas |Tom79|. The optimal interior estimates were obtained by Sogge |Sog: 
3og93]. The same estimate extends uniformly to concave portions of the boundary, 
as shown in |Gri92| for n = 2 and in |SS95 for all n, but not in general (e.g. for 



the 'whispering gallery' eigenfunctions on the disk, see [Gri92] and the remark at 



the end of Section 2.3). The problem of finding optimal L^-bounds for general 



boundary geometry is still open. 

Corollary ^ also follows from the 'sharp Weyl formula' improving 

#{J- A, < A} = 7A/A" + 0(A"-i). 



While this follows from Ivrii's results again, there are simpler proofs, see |See80 
PhaSlj . 



Contents of the paper. In Section ^ we collect some basic facts about our problem 
which are well-known to the experts but scattered or not present in the literature. 
In particular, we prove Corollary |^ from Theorem [l] and give two simple proofs of 
the interior estimate in the case of flat domains. Also, we discuss the torus and 
the disk. In Section we prove estimate outside a boundary layer of width 
A~^, and in Section we derive from this the estimate on u inside this layer, for 
Dirichlet boundary condition. Finally, in Section ^ we describe the modifications 
needed for the Neumann problem. 



2. Basic facts about ||m||oo 

2.1. Multiplicities. The following proposition shows that Corollary ^ is a conse- 
quence of Theorem 0. 

Proposition 3. IfVcz LF'{M) is a subspace of dimension m, then 

(15) sup ||u||oo > |M|-i/2mi/2 

uev 

lhl|2 = l 

where \M\ denotes the volume of M . 

Proof. Let wi, . . . , Vm be an orthonormal basis of V. For simplicity (and sufficient 
for our purpose), we assume that the Vi are continuous, the general case is only 
slightly harder. Define for x,y G M 

m m 

'^yi^) = ^Vi{y)v^{x), a{y) = Uy{y) = ^ |wi(y)P- 

i=l i=l 

We have /^^ a{y) = Y.T=i hi\\l = m, so 

a{y) > m/\M\ 
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for some y. Now by orthonormality of the Vi 

rn 



\U. "2 



y\\2 

1=1 



u 



and \\uy\\oa > Uy{y) = a{y), so HSo/Huy ||5 > a{y) > m/\M\. Therefore 
'I'y/ll^ii/lb satisfies the desired bound. □ 

2.2. The torus. Besides proving Corollary ^ from Theorem |l|, Proposition || has 
another interesting consequence: Let M — R"/(27rZ)" be the 'square' torus. The 
'standard' normalized eigenfunctions are 

(16) Uaix) = (27r)-"/2 gxp(j„ .2,)^ „ g 

with eigenvalue = |ap = af + . . . + a^. These are uniformly bounded, which 
is much better than the bound (^). However, since the multiplicity of is not 
uniformly bounded as A — *■ cxd, one may construct another sequence of eigenfunc- 
tions (as in the proof of Proposition ^ with non-uniformly bounded maxima. In 
fact, the multiplicity of A^ equals the number of representations of A^ as the sum 



of n squares of integers. From standard results on these numbers (see [Gro85 , for 
example) one obtains: 

If M is the square n-torus, then for any A'^ there is a sequence of L^- 
normalized eigenfunctions u with eigenvalues A^ tending to infinity and 
satisfying 

ll^illoo >C^A("-2)/2(log^)iV 

for some Cat > 0. 

As a simple example we take A^ = 5', Z = 1, 2, 3, . . . . This has 4(Z + 1) representa- 
tions as sum of two squares, for example 5 arises from (±1, ±2) and (±2, ±1). This 
gives a sequence as desired for n = 2 and = 1. 



Note also that in the case of the torus one has equality in (113), as follows im 



mediately from ( |16[) . The number-theoretic results referred to above therefore also 
yield the upper bound 

lklloo<aA("-2)/2+^ 



for any e > 0. For more on the case of the torus see Bourgain's article | Bou93 |. Of 
course, the same results hold for a square (or cube) in M" . 



See [Iak97| for more about eigenfunctions on tori and |TZOO for interesting 
recent work on manifolds with uniformly bounded eigenfunctions. 

2.3. The ball. For completeness and to show that the bound in Theorem |^ cannot 
be improved in general, we shortly discuss the case M — {x E M" : \x\ < 1}. In 
polar coordinates, the Euclidean Laplacian is 

(17) A - — + ^^—^— + —As 

where A5 is the Laplacian on the unit sphere S — {\x\ — 1}. By separation of 
variables one obtains that there is a basis of eigenfunctions of the form 

(18) u{ruj) = r-^''-^^/^J„,{\r)<^{Lu), oj e S, 

where $ is an eigenfunction of — A5 (with eigenvalue /i^), J™ is the Bessel function 



(see [VVat48|) of order m = -^//i^ + (n — 2)^/4 (which is always an integer), and A 
is a positive zero of Jm (for Dirichlet boundary conditions). 
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Consider a radial eigenfunction, i.e. — and $ = 1. This is not normalized as 
in (||), so instead of ||w||oo we need to estimate ||u||oo/||w||2. We havej^m ~ (n — 2)/2, 
\\u\\ac = u(0) = C„A^"~^^/^, and the asymptotics 

Jm{s) ^ cs^^l"^ cos(s - rn7r/2 - 7r/4) + 0(s"^/^), s ^ co 
easily imply « /Q\r-("-2)/2(^^)-i/2)2^«-i ^ ^-i^ ||ii||oo/h||2 « A^"-!)/^^ 
and this shows that the bound (g) is saturated by u. 

Similar (but more involved) Bessel function estimates may be used to prove 
Theorem 1^ directly for the ball, and even the stronger estimate for U[a-i,a]- See 



■ Gri92 |. 



Remark: The radial functions are one extreme case of (^^. The other extreme 
case is obtained by taking A = Ami, the first positive zero of J™. For m = 1, 2, 3, . . . 
this yields the sequence of 'whispering gallery eigenfunctions' (say n = 2 for sim- 
plicity). They concentrate on a strip of width k, X^"^!^ at the boundary, as follows 
from (and is made precise by) the estimates 

Ami = TO + am}^^ + 0[vrr^l^)^ as m ^ oo, with a > 

Jmijn + trn^l'^) > m~^/^ for t G [-a/2, a/2] 

Jm{m + tm^/^) < CTO-i/^e^^l*!''' for - to^/^ <t<2a 

for positive constants c, C. (These are easy consequences of well-known asymptotic 



formulas for J^, see [01v54] for example, or |Wat48|, Sec. 8. 4 for weaker but suffi- 
cient bounds.) Sogge showed that away from the boundary concentration can hap- 
pen only on sets of area > A~^/^. This follows from his estimate ||u||6 < CA^/^||u||2. 
In contrast, the whispering gallery eigenfunctions have ||w||6 ~ A-'^/'^||w||2. Note that 
the L°° estimate (|ll) only implies area > A~^ for concentration. This shows two 
things: 

1. Optimal bounds on concentration phenomena are obtained from certain LP, 
p < oo, rather than L°° bounds on eigenfunctions. 

2. As opposed to bounds, these L^-bounds are sensitive to the presence (and 
geometry) of a boundary. 

In general, such optimal bounds are still unknown. 

2.4. General domains in E". Here we prove the interior estimate ( [l^ ) for Eu- 
clidean domains M C M". We give two proofs: First, we finish the wave equation 
proof outlined in the Introduction, and second, we give a more direct proof using 
averaging and Bessel functions. 

First proof: As argued in the introduction, it is sufficient to prove (pj]), with K 
replaced by the wave kernel in R". One way to represent Kjgin is as an oscillatory 
integral, obtained from using the x-space Fourier transform and solving the ordinary 
differential equation that results from (^, 

(19) i^K. {t, X, y) = (2^)-" / cos(ilCI) e^("-^)« d^. 



^ That r "'JmC*") attains its maximum at r = follows immediately from Poisson's integral 
see [Wat48 , Section 3.3. See also Wat48{ , Section 15.31, for a different proof. 



From this it is not hard to get directly. Let us describe a more conceptual 
approach which yields the power n — 1 by pure homogeneity arguments (this is the 



scaling technique in |Mel84| for this special case): From the homogeneity of the 
equation satisfied by Xr™ (or from ( |T^ ) directly) one has the homogeneity 

i^R" {t, X, y) = e^'^Ks.r, (t/ e, x/ e, y/e) 

for e > 0. Since A has constant coefficients, one has furthermore translation in- 
variance K-Rn(t,x,y) = i^n" (t, a; — y, 0). Therefore, (t, a;, x) = K-Rn{t, 0,0) is 
a distribution in t, homogeneous of degree —n. Then its singular support must 
be contained in {0}, and its inverse Fo urier tr ansform is also smooth outside zero 



and homogeneous of degree n — 1 (see | Hor83 |, Theorem 7.1.18), so it satisfies the 



desired bound. Then it is straightforward to see that p * i^Rn " satisfies the same 
bound, i.e. (|Tl|). 

Second proof: Let xo € M, and assume that the ball B of radius R around xq is 
contained in M. We show that 

(20) \uixo)\ < C\^^-'^/-'R-'^^u\\l2^s) 

which clearly implies (p^'). To simplify notation, assume Xq — 0. In this proof, || ||p 
denotes the norm on B, not M. 
Define the spherical average 

considered as function on B. Since h is the average over the functions g*u over 
all rotations g G SO{ti) and since these rotations induce isometrics on L^{B), 
Minkowski's inequality gives 

m2<\\nh. 

Furthermore, h solves (A + A^)/i = (since A is rotation invariant) and is radial, so 
it is of the form (|l8|) with $ = const. Therefore, the same calculation as in Section 



2.3| shows that 

ll^lloc ^ 7->-l/2x(n-l)/2 

Finally, \\h\\oo = \h{0) \ = |u(0)|, so (po|) follows. We remark that this proof may be 



adapted to prove ( [12D instead of (^), see |Gri92|. 

3. Estimates outside the boundary layer 
In this section we prove Theorem |^ for points x with 

dist {x,dM) >X-\ 

As explained in the Introduction, we only have to prove (O) for these x. Propaga- 



tion of singularities (see [EI6r85|) implies that the singular support of the distribu- 
tion t — > K(t, X, x) is contained in the set of lengths of geodesies, possibly reflected 
at the boundary, which start and end at x. Clearly, for small e the only singularities 
in \t\ < e are therefore at t = and possibly at t = ±2dist {x, dM). Therefore, K 
may be expected to be and indeed is representable, for small \t\, as the sum of two 
distributions, a 'direct' term which is only singular at i = 0, and a 'reflected' term. 
To describe the direct term, choose a closed manifold M extending M, and let 

K^"{t,x,y) 

be the solution of the problem (|^) on M. 



To describe the reflected term, it is convenient to introduce geodesic coordinates 
with respect to the boundary; that is, we identify points x of M close to the 
boundary with pairs (x', Xn) € dM x [0, c), c > 0, via the map 

(21) dMx[0,c)^M 

sending {x' , Xn ) to the endpoint of the geodesic of length x„ which starts at x' G dM 
perpendicular to dM. This map is a diffeomorphism onto its image for sufficiently 
small c. We have 

Xn = dist {x, dM). 

Also, let x", a e C, be the distribution on M obtained by analytic continuation in 
a from {Re a > —1}, where it is defined by 



s"/r(a + l) ifs>0 
if s < 



and satisfies {x+^^Y = X+! see [ H6r83 |, Section 3.2. Alternatively, it is the inverse 
Fourier transform of 

-i(a+l)TT/2 }_ 



(22) xU^)=e 

Theorem 17.5.9 in ||Hor85| then gives the following description of the singularity 
of if at i = ±2dist {x,dM): 

Proposition 4. For sufficiently small e, there is a distribution 

I{x',e,t) e V'{dM X E X (-e,e)) 
so that for \t\ < e and x„ > we have 

2t 

K{t, X, x) = X, x) - t~"I{x\ -Y',t). 

Furthermore, I has .support in \6\ < 1 and .singular support in \6\ = 1, and near 
9 = 1 we have 

N-l 

I{x',9,t) = a,(a:',^?,i)x+"^"^'^^'(l - 0) + Rn{x' ,e,t) 
j=o 

with smooth functions aj and a continuous remainder R^, for any N > (n — l)/2. 
A similar expansion exists near — —1. 



Note: In the statement of Theorem 17.5.9 in |H6r85 the direct term (called 
t~'^Ii there) is not identified like in the statement above. But in the proof (middle 
of page 59, loc. cit.) it is chosen as K'^" like above. 

The proposition is obtained by analyzing the regions t <3a;„ and t > 3x„ sepa- 
rately: In t < 3xn one uses a suitable scaling (cf. Section 2.4) and the Hadamard 
parametrix, and in t > 3x„ (where no singularities should occur) one uses propa- 
gation of singularities estimates (this is the harder part). 

Since the estimates for K'^" are already known by the interior estimate, (|l^) will 
follow from: 

Lemma 5. For Xn > A^^, we have 

I / e''^p{t)t-"I{x',—,t)dt\<CX''-\ 
Jo t 



For convenience, we restrict to positive t. Negative t are handled in the same 
way. 

Proof. The integrand is zero for t < 2a;„ and for t > e. We spht the integral up 
into a part where t > 3a;„ and a part where t < 4a;„ , using a cutoff function smooth 
in t/xn- By Proposition^, / is smooth on the first part and therefore bounded, so 
this part is dominated by a constant times 

t^"dt<x-"+^ < A"-^ 

In the second part, where t < Axm we split up / as in the proposition, with a fixed 
N > {n — l)/2. Then the term i?jv can be handled in the same way as the first 
part. 

It remains to analyze the singular terms. If we denote the cutoff function by 
ip{t/xn), tp G C'^(l,4), they take the form (assuming for simplicity that p was 
chosen to be constant near zero) 

e^^\{x',^,t)xl(l-^m-)t-''dt, 

a = j — {n + l)/2. If we change variables r = t/2xn and use the homogeneity of 
X+, this becomes 



x',T,t)xl{T-l)dT 



with b smooth and supported in r S (1/2,2). Using ( |22| ) one sees by a short 
standard calculation that this is bounded by a constant times 

a;-"+i(Ax„)("-i)/2-^' = A"-i(Ax„)-("-i)/2-j < 

□ 

4. Estimates in the boundary layer 

So far, we have proved \u{x)\ < CA^"-!)/^ for dist {x,dM) > A'^ The proof of 
Theorem |l| in the Dirichlet case will therefore be completed by the following lemma: 

Lemma 6. If u is a solution of 

{A + X^)u = 
vanishing at the boundary of M , then 

(23) max |u(a;)| < max |w(x)|. 

x:dist(x,dM)<\-'^ x:dist (x,dM) = X-^ 

Proof. Without loss of generality, assume u is real valued. Under the identification 
near the boundary given by (^l|), the metric takes the form 

g{x) = dxl +g'{x',Xn) 

where g'{-, Xn) is a Riemannian metric on dM for each Therefore, the Laplacian 
has the form 

(24) A = — + a{x)—+Pix',Xr,,D,,) 

C/t/^yj KJ J1 
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with a Laplacian P on the boundary (depending on the parameter a;„) and a smooth 
function a. The generahzed maximum principle says that if w is a positive function 
on the strip 

S ^{0<Xn< A"i} C M 

with 

(A + A2)w < 

then 

\u\ \u\ 

(25) max— < max — . 

S V OS V 



(Apply |PW84], Theorem 10, to u and then to —u.) We apply this with the function 



TT 3 

v{x',Xn) = sin(- + -{Xxn - !))• 
We have v > 0.07 on S and thus, from (p4|), 

(A + X^)v = -^X^v + 0(A) < on 5 for large A, 
so the generalized maximum principle applies. Now m = at the outer boundary 



a;„ = of S, and u = 1 at the inner boundary x„ = A ^. Since w < 1 on S*, (25) 



implies (|2^). □ 

5. The Neumann problem 

Theorem |^ and Corollary ^ remain true if the Dirichlet boundary condition 
u^gj^j = is replaced by the Neumann boundary condition 

dnUlOM = 

in (|l|), where 9„ denotes the outward normal derivative. 

Let us sketch the proof: The proof of Proposition ^ (and thus the proof of 
Corollary |^ from Theorem [^) as well as the reduction of Theorem ^ to the wave 
kernel estimate ( pj| ) carry over literally, except that K is replaced by the Neumann 
wave kernel (i.e. dnK^^^QM = instead of K^qm = in (^, where dn refers to the 



x-coordinates). (Also, Section 2.4 applies literally to the Neumann problem.) 
Proposition^ holds for the Neumann wave kernel as well (with different aj), by 



straightforward modification of the arguments in |Hor85 | 



Finally, Lemma |^ remains true for Neumann eigenfunctions, except that ( |23| ) 
must be replaced by 

(pSf) max |m(x)|<20 max |u(a;)|. 

' — ' x:dist{x,dM)<X-^ x:dist(x,dM)=X-^ 

To see this, just note that, in the comparison of the functions u and u, the maximum 
of u/v over S not only must occur at a boundary point xq of S', but also the outward 
normal derivative at xq must be strictly positive. See Theorem 10 of | PW84| | (this 



is sometimes called Zaremba's principle). Now choose v — sin(^ + |Ax„), then 

dnV — ~dv/dxn ~ for Xn = 

and so dn{u/v) = for Xn = 0. This means that the maximum of u/v over S must 
occur at the inner boundary a;„ = A~^ of S. Since < 20 there, we get (pBf) 
after applying the same argument to —u. 
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